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Conserved quantities are crucial in quantum physics. Here we discuss a general scenario of Hamil¬ 
tonians. All the Hamiltonians within this scenario share a common conserved quantity form. For 
unitary parametrization processes, the characteristic operator of this scenario is analytically pro¬ 
vided, as well as the corresponding quantum Fisher information (QFI). As the application of this 
scenario, we focus on two classes of Hamiltonians: su(2) category and canonical category. Several 
specific physical systems in these two categories are discussed in detail. Besides, we also calculate 
an alternative form of QFI in this scenario. 
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I. INTRODUCTION 

With the development of quantum physics, quantum 
mechanics nowadays is way beyond the phase of logical 
argument and mind experiments, but deep into the ap¬ 
plied field and even our daily lives. Quantum metrology 
is one of the most successful applications of quantum me¬ 
chanics. Starting from the pioneer work of Caves mm, 
people gradually realized that many quantum effects and 
states are available to enhance the measure preci¬ 

sion of physical parameters via various quantum systems. 

Quantum Fisher information (QFI) is a central quan¬ 
tity in quantum metrology because it depicts the lower 
bound on the variance of the estimator 0 due to the 
Cramer-Rao theory: S 2 0 > (zaF) - 1 [161, UJ , where 5 2 0 
is the variance, v is the number of repeated experi¬ 
ments, and F is the QFI. Generally, the QFI is defined 
as F = Tr (pL 2 ). Here p is a parametrized state and L 
is the symmetric logarithmic derivative (SLD), which is 
determined by the equation dop = (pL + Lp)/ 2. 0 is 
the parameter under estimation. Quantum Fisher infor¬ 
mation matrix (QFIM) T is the counterpart of QFI in 
multi-parameter estimation. The element of QFIM is de¬ 
fined as Fxxin — Tr(p, L n j-)/2. T m , L n are the SLD 
operators for rath and nth parameters under estimation, 
respectively. 

Because of the importance of QFI, its calculation is 
always an interesting and attractive topic. Recently, 
several alternative formulas of QFI have been devel¬ 
oped [13 HO HMZ1- For unitary parametrization pro¬ 
cesses, new expressions of QFI and QFIM were given for 
both pure and mixed states [201 HD- The QFI and QFIM 
here are totally determined by a characteristic operator 
TL and the initial states. In these expressions, all the 
information of parametrization is involved in TL. Fur¬ 
thermore, TL can be written into an expanded form m, 
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which is particularly useful when the nth order commu¬ 
tation between the Hamiltonian and its partial derivative 
can be truncated or is periodic. Moreover, for a general 
parametrized exponential state, the SLD operator is also 
found to be expressed in an expanded form [23]. These 
facts prompt us to study various Hamiltonians owning 
above properties. This is the major motivation of this 
paper. 

Conserved quantities are crucial in both classical and 
quantum physics. Since Emmy Noether connected them 
with differentiable symmetries in Noether’s theorem, 
searching for conserved quantities becomes a prior mis¬ 
sion when facing a novel new system. Locating conserved 
quantities will not only help us to find hidden symmetries 
of these systems, but also give us a easy perspective to 
describe and classify them. 

In this paper, we discuss a scenario, in which all sys¬ 
tems share a common conserved quantity form. For the 
unitary parametrization processes, we provide the ana¬ 
lytical expressions of the characteristic operator and the 
QFI. The maximum QFI and the corresponding opti¬ 
mal initial states are discussed. Moreover, we also study 
QFI for the parametrized thermal states in these sys¬ 
tems. The scenario we discuss includes many systems, 
of which two typical classes, su(2) category and canon¬ 
ical category, are the main applications of this paper. 
In the su(2) category, ferromagnetic two-spin system, 
anisotropic two-spin system and a spin-one system are 
discussed in detail. In the canonical category, a cavity 
optomechanical system are provided as an example. At 
the end of this paper, we discuss an alternative form of 
QFI given by Luo et al. [28j. The formula of this alter¬ 
native QFI for above scenario is analytically calculated 
and discussed. 

The paper is organized as follows. In Sec. [TT| we re¬ 
view quantum metrology with unitary parametrization 
processes. In Sec. |III[ we propose a scenario in which all 
systems share a common conserved quantity form. For 
unitary parametrization processes, the characteristic op¬ 
erator and the QFI are analytically provided. We also 



2 


discussed the QFI for the parametrized thermal states 
in the scenario. In Sec. the applications of this sce¬ 
nario, including two typical classes: su(2) category and 
canonical category, are given and discussed. In Sec. [V| 
an alternative form of QFI is calculated for this scenario. 
Section |VT] is the conclusion of this paper. 


II. UNITARY PARAMETRIZATION 

The unitary parametrization process is a widely used 
parametrization strategy in quantum metrology. Re¬ 
cently, it has been found that the QFI and QFIM for 
both pure and mixed states can be expressed via a char¬ 
acteristic function H I2DH2H- Denoting the parameter 
under estimation as 6, H is defined as mm 

He ■■= i (deU^) U, (1) 

where V is a unitary parametrization transformation, i.e., 
pe = UpoW with po a 0 -independent density matrix. 
Now we consider the situation that the transformation is 
generated by a time-independent Hamiltonian, namely, 
U can be written in the form 

U = exp [—itH{6)\, (2) 

where 6 = (#i,# 2 >---) T is a vector of parameters under 
estimation and the parametrized Hamiltonian H(6) is 
time-independent. Based on a recent work [21], the char¬ 
acteristic operator He for parameter 0 can be expressed 
in an expanded form 


where (A 2 'Hg) i := (H 2 e )i - (He) 2 is the variance of the 
characteristic operator on |^). For a purely initial state, 
the element of QFIM reduces to m 

mn = 4COV (Hmi Hn) • (7) 

The covariance is taken on the initial state. Based on 
this equation, the QFI for purely initial states is [20, 21 j 

F e = 4(A 2 H e ), (8) 

namely, the QFI for parameter 6 is actually the variance 
of the corresponding characteristic operator. 

In the expression of QFIM in Eq. for any m, if 
we assume H m = H' rn + c m , with c rn a complex number, 
it can be checked that co\ i(H m ,H n ) = covi(H f m , H' n ). 
Meanwhile, in the second part of Eq. the overlap 
(^i\cm\^j) always vanishes for any m when i ^ j. Thus, 
Hm and H ^ share the same expression of QFIM. This 
fact indicates that the number terms of characteristic op¬ 
erator do not affect the value of QFIM, therefore they can 
be neglected in metrological problems. 


III. CONSERVED QUANTITIES 

Conserved quantities are important in theoretical 
physics. In the following we propose a scenario in which 
all systems share a common conserved quantity. At first, 
we introduce a ^-dependent Hermitian operator V, which 
is defined as 


(j+\ n + 1 

’ H ‘ = i E(rYy.^) nd > H ’ ( 3 ) 

n=0 ' ' 

where H x = [H, •] is a superoperator. With a known H, 
the QFIM can be obtained easily. The element of QFIM 
can be expressed by [21] 

M 

mn ^ ^ IpjCOV^ ( Hmi'H'n ) 

_ £ ^A Re MilHrnWjWjWnWi)), (4) 
i& Pi+P 3 

where the covariance reads 

COV^ (Hm,Hn) = ('Hm)i^n)i- (5) 

In above equations, H m is short for Ho m . Pi and |^) 
are the ith eigenvalue and eigenstate of initial state po, 
which means pi and | are independent of 0. M is the 
dimension of the support of po. (•)$ is the expected value 
on |^), i.e., (-)i = ( / ijji | • | ipi). It is known that the diag¬ 
onal elements of QFIM are the QFIs for corresponding 
single-parameter estimations, thus, the QFI for a unitary 
parametrization process can be written into the form [21] 

m 8 

Fe = J2±Pi(^ 2 Kg) i - F (6) 

<=1 -y Pi ' Pj 


V = 


(h x ) 2 -n 2 


d e H, 


(9) 


where U 2 is a real number and 0 is a parameter in the 
Hamiltonian. This operator is generated by the Hamil¬ 
tonian H and related to the parameter 6. The central 
quality of the scenario we discuss here is that V is a con¬ 
served quantity for all systems in this scenario, namely, 
the Hamiltonian H satisfies 


[V,H}= 0. (10) 

Based on the definition of V, above equation can be 
rewritten into 



H x d e H = 0. 


(ii) 


This equation implies that H x 8qH is the eigenopera- 
tor of superoperator ( H x ) 2 1 with U 2 the corresponding 
eigenvalue. From this aspect, one can check that H x 8qH 
and ( H x ) 2 8qH are also eigenoperators of (H x ) 2n with 
U 2n the eigenvalues for n > 0, i.e., 

[(H x ) 2n - n 2n ] (H x )%H = 0, ( 12 ) 

where z = 1,2. In the following we take 6 as the param¬ 
eter under estimation. Based on this equation, the char¬ 
acteristic operator He in Eq. © will be separated into 
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two parts via the parity of n. Through some straightfor¬ 
ward calculations, the analytical expression of He can be 
obtained as below 

Ue = [-t-i(d t f)H x + f(H x ) 2 ]d e H. (13) 

Here / is short for the function /(U,£), which is defined 
as 


/(ft, t) ~ [ftf - sin (Hi)]. (14) 

If D > 0, / is positive and monotone increasing with the 
passage of time. Utilizing the conserved quantity V, the 
expression of He can be rewritten into 

H e = fV + (/ft 2 - t) d e H - i ( d t f) H x d e H. (15) 

Moreover, from the expression of /, the coefficients can 
be simplified as /U 2 — t = —U _1 sin (fit) and 8 t f = 
2fi -2 sin 2 (fit/2), then F# can be finally expressed by 

He = fV- / sin(fti)3*.ff - ^ sin 2 (^i) H * d e H. (16) 


This is the general expression of the characteristic op¬ 
erator He for the scenario in which V is a conserved quan¬ 
tity. In this scenario, the characteristic operator is the 
linear combination of the conserved quantity V, partial 
derivative 3qH and commutation between 8qH and H. 
With the expression of He, the QFIM and QFI can be 
calculated through Eqs. @-(§ for both mixed and pure 
states. For the long-time limit or the situations that the 
value of U is very large, the oscillating terms in Eq. (16) 
can be neglected and He reduces to 


= ft* V. 


(17) 


The characteristic operator is then proportional to V. For 
purely initial states, the QFI is 


F e 


^< A2 V), 


(18) 


where the equality {d e H,H x d e H} = H x (d e H) 2 has 
been applied. {■, •} represents the anticommutation. A 
simple example of this scenario is that H x 8qH is pro¬ 
portional to 8qH, i.e., H x d$H = TtdeH. In this example, 
V = 0, then H reduces to iQ~ 1 (e znt — 1 )8qH [2T]. Espe¬ 
cially, if U can be chosen as iuo with u a positive number 
in this example, the characteristic operator He will re¬ 
duce to —deH/ue for the long-time limit. A realistic case 
of this scenario is a collective spin system in an external 
magnetic field, which has been discussed in Ref. m in 
detail. 

Moreover, if the operator Vi = ( H x — f T)8eH is a non¬ 
trivial conserved quantity, then V = DVi is also a non¬ 
trivial conserved quantity. Therefore the characteristic 
operator in this case can also be expressed in the form of 
Eq. Jl6). 


Thermal states of the scenario .-Thermal states widely 
appear in realistic world. In quantum theory, a general 
thermal state can be written as 


p = exp {-PH ), 


( 21 ) 


where /? = l/{k&T) and the partition function Z = 
Tv{e~ 0H ). T is the temperature and k b is the Boltzmann 
constant. In Plank unit, &b = 1. Recently, Jiang |23| 
provides the expression of SLD for a general exponential 
state pe = exp[G(0)], which is 




v2 n 


n =0 


(2n + 2)! 


8 0 G. 


( 22 ) 


Here £>2 n +2 (2 n + 2)th Bernoulli number. For 

the thermal state expressed in (21), G = —/3H — In Z. 
Then it is easy to check the equality (G x ) 2n 8eG = 
(— P) 2n+1 (H x ) 2n 8eH . For any system in the scenario, 
V is a conserved quantity. The SLD operator can then 
be calculated as 


L = p[rV+ (rft 2 - 1) d e H] , 
where r is short for r(/3, fi), which is defined as 


(23) 


namely, it is actually determined by the fluctuation of 
the conserved quantity. 

When V is a number (which is a trivial conserved quan¬ 
tity), it can be neglected according to the analysis in 
Sec. [IIJ Thus, the characteristic operator in this case re¬ 
duces to 


r{P,n) 


1 

U2 


1 _ /3ft tanh 


(24) 


The regime of r is [0,ft -2 ]. Meanwhile, one can see that 
rft 2 — 1 = — tanh(/?ft) /(/3ft). Then the SLD can be al- 
ternatively written as 




i2 


sin (ftf) d e H - — sin 2 ( -t )H x d e H 


ft 


• ( 19 ) 


For purely initial states, the corresponding QFI is 


Fe ft 2 L 


sin 2 (nt) ((deH) 2 ) - -^sm 4 l^t)((H x d e H) 2 ) 


i 4 


ft 


-—sin(ftf)sin 2 ( —t ) (H x (d e H) 2 ) 


( 20 ) 


L = /3rV- A tanh (/3ft) d e H. (25) 

With above equation, the QFI Ft = (L 2 ) can be finally 
expressed by 

F t = /3 2 r 2 (V 2 ) T + T tan h 2 (/3ft) ((d e H) 2 ) T 
— ft r tanh (^) ({V,d e H}) T . 


(26) 
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Here (-)t is the expected value on the thermal states. 
If the thermal states in Eq. (21) can be rewritten into 
the form Ue~^ H °U^ with Hq a parameter independent 
Hamiltonian, the QFI can also be calculated utilizing the 
function He in Eq. (16). At the zero-temperature limit, 
tanh(/?fi) = 1 and r = D -2 , the SLD reduces to L = 
/3Fl~ 2 V and the QFI is F T = /3 2 Ft~ 4 (V 2 ) t- Similarly, at 
the high-temperature limit, tanh(/3U) ~ /3Q and r = 0, 
the SLD is L = —(3doH and the QFI can be written as 
F T = /3 2 ((d e H) 2 ) T . 


IV. APPLICATIONS 



In the following we will solve the metrological prob¬ 
lems in some realistic systems in the scenario where V 
is a conserved quantity. We mainly focus on two classes 
of Hamiltonians. The first one is mainly related to the 
generators of su(2) algebra and is called su(2) category; 
the second one is related to the canonical variables: the 
position operator x and the momentum operator p, and 
it is called canonical category. We first discuss su(2) cat¬ 
egory. 


A. su(2) category 


Figure 1: (Color online) Optimal points to access the maxi¬ 
mum quantum Fisher information for ferromagnetic two-spin 
system at the long-time limit. The solid blue, dashed red, 
dash-dot black and dotted pink lines represent the optimal 
points for B = 0.4, 0.8, 1.0 and 10, respectively. 


Utilizing above two equations, we can obtain 



4 (1+4 B 2 f 


H?d e H 1 = 0 . 


(31) 


This equation implies that if we choose 


Ferromagnetic two-spin system 


fi = 2-/L + 4B 2 , (32) 


As the first application of the su(2) category, we now 
consider a ferromagnetic two-spin system in an external 
magnetic field. The Hamiltonian of this system reads 

Hi = -afaf-B(<7f + <7|), (27) 

where a\ = cq 0 11 and o\ = 1 0 o i for i = x, y, z. 
(Ji is a Pauli matrix, t is the identity matrix, and B 
is the strength of the external field. The optimization 
of QFI in a general Ising model with GHZ-type state 
has been discussed recently [22]. In this case, we take 
B as the parameter under estimation. Before the main 
calculation, we introduce three operators 

J X = \ (o-fo-f - of 0 !) > 

< Jy= I (of O'! + < 7 M)> ( 28 ) 

Jz= iK + ^)- 

It is worth to notice that J x , J y and J z satisfy the su(2) 
commutation [Ji,Jj\ = FijkJk with the Levi-Civita 
symbol. In addition, the anticommutation is {Ji,Jj} = 
2 5ijJi with Sij the Kronecker delta function. Using these 
operators, the first and second order commutations be¬ 
tween the Hamiltonian and its derivative are calculated 
as below 


H?d B H x = -i8J y , (29) 

{H?fd B H i = 16(2 BJ X -J Z ). (30) 


the operator V defined in Eq. © is a conserved quantity. 
Specifically, it is 


V = 32 Bv-J, 


(33) 


where J = (J x , J y , J Z ) T and v = (1,0,2 B) T . Using this 
method, we find a non-trivial conserved quantity in this 
two-spin system. Moreover, based on the property of 
quantum conserved quantity, all the operators, for which 
the corresponding vectors share the same or opposite di¬ 
rections with u, are conserved quantities, v is a vector in 
the x — z plane. When B is zero, v is along the x axis. 
With the increase of H, v rotates around the y axis from 
x axis to z axis. For a very large L>, v is almost along 
the z axis. 

Compared with the expressions of Q and u, one can 
see that Ft = 2v with v = \v\. the function / can then 
be rewritten into / = [2 vt — sin (2 vt)\ /8v 3 . The charac¬ 
teristic operator Hb can be calculated via Eq. (13). Its 
explicit expression is 


(34) 

where the vector x reads x = (8L>/, — 2d t f, t — 4/) T . With 
above Hb , the QFI for purely initial states in this case 
can be written as 


F b 


16 

y 


|f| 2 (|jT)-3(f-(J)) 2 


(35) 
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where |J| 2 = 3(1 +erf erf )/8. During the calculation of 
Fg, the relation {Ji,Jj} = 2 SijJi has been used. Above 
equation of QFI implies that the its maximum is attained 
when (| J| 2 ) is maximum and (J) is vertical to +. In a 4- 
dimensional Hilbert space, the maximum value of (erferf) 
is 1, which indicates the formula of the maximum QFI 
must be 


F Bi max = 4|^| 2 , (36) 

and the optimal initial state is required to be in the form 

|V’o P t)=oi|00)+o 2 e^|ll), (37) 

where is a real number. To make x • (J) = 0, the 
amplitudes need to satisfy the equation 

CL 1 CL 2 (x x cos p + x y sin p) + -x z (a\ — a|) = 0. (38) 

Here Xi (i = x,y, z) is a element of x. Since x z = t — 4/ 
is not always zero for t > 0, then <21 <22 = 0 cannot be 
a solution of this equation, which means above equation 
can be further simplified into 


x x cos p + x v sin p + \x z 




(39) 


All states satisfying this equation are available to access 
F^max- At the long-time limit, / ~ t/Av 2 and d t f <C /, 
then Eq. (39) reduces to 


4 B 


cos p + B 




(40) 


Figure [l] gives the optimal points (0, ai/a 2 ) to access 
the maximum QFI for different B at the long-time limit. 
The solid blue, dashed red, dash-dot black and dotted 
pink lines in this figure represent the optimal points for 
B = 0.4, 0.8, 1.0 and 10, respectively. From this figure, 
it can be found that with the increase of F, the curve 
becomes more flat. This behavior indicates that for a 
strong external field, the maximum QFI is insensitive to 
the relative phase <p. This is actually due to the fact 
that when B is very large in above equation, cosp can 
be neglected and the equation reduces to 


<21 a 2 
a 2 a\ 


(41) 


One solution of this equation is a\ = a 2 . Thus, for a 
strong external field, the maximum QFI can be saturated 
at <21 = <22 with any phase p. 

A widely studied special form is that a\ = = l/\/2, 

i.e., IV^opt) = (|00) + e^|ll))/V% then the optimal rela¬ 
tive phase p opt reads 



Figure 2: (Color online) The optimal relative phase p opt as a 
function of B and t for the optimal state (|00) + e l( ^\ll))/y/2. 


field. This fact means at the long-time limit, the state 
(|00) + e z ^|ll))/\/2 is always an optimal state for any 
strength of external field. Figure [2] shows the general 
variation of optimal relative phase as a function of B 
and t. For a weak external field, p opt is growing rapidly 
with the increase of B and t. This is because when 
the external field is very weak, 4 Bf/d t f — 4:Bt/3 , then 
p opt arctan(4F£/3) ~ 4F£/3. Thus, the optimal phase 
Popt grows almost linearly with B and t in this regime. 
With the continue increase of B and £, p opt shows a os¬ 
cillating behavior for intermediate strength of external 
field. The oscillation amplitude of p opt trends to shrink 
with the passage of time as p opt = 7t/2 when t is infinite 
large. 

Thermal state- Here we consider the thermal state of 
this ferromagnetic two-spin system, of which the den¬ 
sity matrix reads pi = exp(— f3Hi)/Zi, with Z\ = 
Tr [exp(— /3Hi)\. It is already known that V is a conserved 
quantity, then based on Eq. (23), the SLD operator for 


B of this thermal state can be written as 


Lb — PrV H— tanh (2 f3v) J z 
v 


(43) 


where r is defined in Eq. (24) and in this case, it has the 
form 


r = 


4u 2 


1 ~ 2 \j3v tanh ^ 2 ^ 


(44) 


Based on the expressions of Fb, the QFI can be calcu¬ 
lated as 


F t = 2 P 2 (16 v 2 r 2 - 8r + l) (1 + (oferf)). (45) 


Popt = arctan 



(42) 


Ft here is determined by the correlation function (erferf), 
which can be analytically solved as 


At the long-time limit, this optimal phase reduces to 
a constant number 7r/2, independent of the external 


(of 0 -f) = -1 + 


2 cosh (v/3) 


cosh (v/3) + cosh /3 


(46) 
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Figure 3: (Color online) Quantum Fisher information as a 
function of B for thermal state of the ferromagnetic two-spin 
system. The temperatures are set as T — 0.5, 1.0 and 1.5 for 
the solid blue, dashed red and dash-dot black lines. 


Thus, Ft can be finally expressed in the form 


F t = (l 6 u 2 r 2 — Sr + l) 


4 /3 2 cosh (v/3) 
cosh (v/3) + cosh (3 ’ 


(47) 


Figure [3] shows the variation of Ft as a function of B. 
The temperatures are set as T = 0.5, 1.0 and 1.5 for the 
solid blue, dashed red and dash-dot black lines. In this 
figure, for the weak external field, the QFI grows greatly 
with the increase of B. However, for the strong external 
field, this growth is not significant. With respective to 
the temperature, the decrease of T shows a positive effect 
on the QFI and is useful for the precision measure in this 
system. 

For a very low temperature and nonzero external field, 
cosh/3/cosh(u/3) ~ ~ 0, then the QFI in Eq. (47) 

reduces to the form 


(48) 

This equations shows that in the low-temperature regime, 
the QFI is inversely proportional to the square of T. 
Thus, the decrease of T will dramatically improve the 
value of Fp. Meanwhile, the increase of B will also en¬ 
hance the value of QFI. The maximum value of Ft in 
this regime is F T>max = 4/T 2 . 


Anisotropic two-spin system 

In the following we consider a more general case: the 
anisotropic two-spin ferromagnetic XY model with an in¬ 
homogeneous external magnetic field. The Hamiltonian 
of this system is 

tt _ 1 A T x x 1 T y y 

2 — ^ a l°2 2 a l a 2 

-B+ (erf + <rf) - B- (af - erf), (49) 


where 7 is the anisotropic parameter. This Hamiltonian 
can reduces to the Hamiltonian Hi in Eq. (27) with 7 = 1 
and F_ = 0. Before the main calculation, we introduce 
a new group of operators 


(s x =\ (a?a$ - a\al ), 

\Sy= 2 (of of + ofof), (50) 

[Sz= \ (o"f of) • 


Similarly with J x ,y,z , S x ,y,z also satisfies su( 2 ) commu¬ 
tation, i.e., = ieijkSk • Meanwhile, the anti¬ 

commutation relation is {Si,Sj} = 2 SijSi. A very in¬ 
teresting property between these two groups of operators 
is that 


JiSj = SjJi = 0, Vi, j = x, 2 /, z. 


(51) 


[Ji,Sj] = 0 for any i and j is a natural result of this 
property. 

Utilizing these two set of operators, Hamiltonian (49) 
can be written as the sum of two parts, i.e., 


F 2 = 2(F + + F_), (52) 


where the sub-Hamiltonians H+ is only related to J x , y , z 
and H _ is only related to S XrVrZ . Their specific formulas 
are 


F+ = - 1 J X -2B + J Z , (53) 

H_ = -S y -2B_S Z . (54) 

Since Ji and Sj are commutative for any i and j, H+ 
and H_ are also commutative. Here we take both F + 
and B _ as the parameters under estimation. Through 
some algebra, we find that if one choose 

= 2^/ 7 2 + 4B 2 , = 2/mF, (55) 

the corresponding operator V + and V_ defined in Eq. © 
are conserved quantities and have the form 

V+ = 32 B + v+ • J, V_ = 32 F_u_ • S, (56) 

where J = ( J x , J y , J Z ) T , S = ( S x , S y , S Z ) T and the vector 
and v _ read 


v+ = ( 1 ,0,2B + f, v- = (0,1, 2£?_) t . (57) 

Then there is Q± = 2v± with v± = |F±|. 

Utilizing these conserved quantities, the characteristic 
operators F+ and H- can be expressed by 

F+ = 4x+ • J, H-= 4x- • S. (58) 

The vector oc± in above equation reads 


x + = (8'iB + f + ,-2'jdtf+,t - 47 2 /+) T , (59) 
x- = (25 t /_,8B_/_,i-4/_) T . (60) 
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Figure 4: (Color online) Optimal points to access the max¬ 
imum quantum Fisher information for anisotropic two-spin 
system at the long-time limit. The solid blue, dashed red, and 
dash-dot black lines represent the optimal points for 7 = 0.3, 
0.6 and 0.9. 


Here f± = [2v±t — sin (2v±t)\ /8v±. For purely initial 
states, the QFI is the variance of TL on the initial state. 
For parameter H + , the QFI F + shares the same form 
with that in Eq. (35), i.e., 


For the long-time limit, this equation reduces to the same 
form with Eq. (40). Moreover, taking b\ = 62 = l/\/2, 
the optimal relative phase can be written as 


^opt = arctan j • ( 66 ) 

At the long-time limit, ip opt also equals to tt/2. 

Now we consider the situation that both and B_ 
are unknown parameters simultaneously. It is known 
the condition for the saturation of the multiparameter 
Cramer-Rao bound in unitary parametrization process 
is [ 21 ] 


(^in|[H + ,H_]|^in)=0. 


(67) 


From the expressions of TL± in Eq. ( [58] ) and the property 
that Ji and Sj are commutative for all i and j, it is easy to 
see that [%+, TL-\ = 0. Thus, and B_ can be jointly 
measured for any purely initial state. Since TL+ and TL- 
are commutative and based on Eq. (51), the off-diagonal 
element of QFIM is 


?+- =?-+ = -m+)(n-) 


( 68 ) 




|f + | 2 (|J| 2 )-3(f + - 


(61) 


For parameter H_, the QFI F_ has the similar form 


F_ = 


16 


|f_| 2 (|5| 2 )-3(f_-(5)) 2 


(62) 


where |Sj 2 = 3(1 — crfcr^/S. The maximum value of 
F + is 4|7f + | 2 and the corresponding optimal initial state 
has the same form as that in Eq. (37). For the long-time 


limit, the equation for B+ that the optimal points satisfy 


is 


4y B+ cos + (ui — 7 2 ) (— — —^ = 0. (63) 

\ a 2 0 * 1 / 

Figure [4] gives the optimal points to access the maximum 
QFI for the long-time limit. The solid blue, dashed red, 
and dash-dot black lines represent the optimal points for 
7 = 0.3, 0.6 and 0.9. It is found in this figure that with 
the increase of 7 , the curve of the optimal points gets 
more sharp, indicating that the maximum QFI is more 
sensitive to the relative phase for a large 7 . 

Similarly, the maximum value of F_ is 4|x _| 2 with the 
optimal initial state 

l^opt) = bi|01) + bae^llO), (64) 


The expected value above is taken on the initial state. 
According to the Cramer-Rao theory, we have 

5 2 B ± > - , (69) 

J- + + ./* ~ J ~+- 

When the initial state is chosen as |Vy>pt) or |4> opt ), 
T+- vanishes and the inequality above reduces to the 
form of the single-parameter cases, which indicates that 
the joint measurement of and B _ can be performed 
by |Vk>pt) or |$ opt ). However, it should be noticed that 
l^opt) and |4> 0 pt) are orthogonal, which means even B+ 
and B_ can be jointly measured, there does not exist an 
optimal state to access the maximum QFI for both 
and B_ simultaneously. This fact implies that the joint 
measurement here is not as good as the single measure¬ 
ment. 

Thermal state .-Next we consider the thermal state of 
# 2 , which is p 2 = exp(—/ 3 i 72 )/^ 2 , with Z 2 the partition 
function. Based on Eq. ([23]), the SLD operators for B+ 
and B _ can be expressed by 


L± = /3r±V± -tanh ( 2 /3v±) d±H±. (70) 

v± 

where d± is short for 8b± and the coefficient 


r± = 


4vl 


1 - 


2 j3v± 


tanh ( 2 /3v±) 


(71) 


where b\ and 62 are real numbers. The equation for 61 , 
62 to satisfy to access the maximum QFI is 

2 d t f- sirup—8B-f- cos<p=^(t- 4 /-)^ “ ( 65 ) 


Since {^, Sj} = {J^Jj} = 0 for i ^ j, the thermal QFI 
for and T>_ can then be written as 

Ft+ = 2/3 2 (l67 2 i; 2 r+-87 2 ?'+ + l)(l + (7 cr 2)) I (72) 
Ft- = 2/3 2 (l 6 v 2 r 2 - 8r_ + l) (1 - <crfcr|>) . (73) 
















-10 0 
0 0 0 

0 0 1 


(81) 


For the thermal state, the correlation function (o-fo-f) is 


(crfo-f) = -1 + 


2 cosh (pv+) 

cosh (/3v+) + cosh (/3v-)' 


In the low-temperature regime, there is 


cosh(/fa>_) ^ e p (v _- v+) 
cosh (/3v+) 


(74) 


(75) 


Therefore, when v+ equals to V-, e^^ v ~ equals to 1, 
and (erf erf) = 0 . Fp,+ and Ft,- reduce to 



When v+ is smaller than r>_, e f3 ( v -~ v+ ' ) trends to infinity 
in the low-temperature regime, then (erf erf) — —1 and 
Ft,+ and Ft,- is in the form 


F t>+ ~ 0, Ft - ^ T . (77 ) 


When v+ is larger than r>_, e^ v - ~ 0, then (erferf) — 
1 , and the QFI Fp,+ and Ft,- can be written as 

F T,+ - f2 ’ Ft ’~ ~ °- ( 78 ) 

Above analysis shows that when v+ is smaller than 
V- 7 the parameter F + can be barely estimated via the 
Cramer-Rao theory, so as F_ when v+ is larger than V- . 
Thus, if either of F + and F_ is the parameter under es¬ 
timation, we have to tune down the value of the other 
one to make sure that nonzero QFI exists. 


Spin-one model 


Not only the spin-half systems, but also some spin-one 
systems, can fit in the su( 2 ) category. Here we show such 
a spin-one system in the one-axis twisting model. The 
Hamiltonian of a one-axis twisting model with a trans¬ 
verse field can be written in the form [[29H32] 

Hs = xJx+BJ zy (79) 

where J x = (a^b + atf)/ 2 and J z = (a) a — b^b)/2 are the 
Schwinger operators. Another one is J y = {a)b — ab^)/2i. 
X is the coupling constant and B is the strength of the 
transverse field. This Hamiltonian can be realized in 
many physical systems including two-component Bose- 
Einstein condensates |33j[34]. Now we consider one real¬ 
ization that a two-boson system in a double well. Since 
the particle number is a conserved quantity, this system 
can be expanded in the basis {|02), 111), |20)} in Fock 
space, which can be mapped as a spin-one system. In 
this basis, the Schwinger operators has the form 


Jx 


1 (° 10 
— 1 0 1 

^\0 1 0 


Jy — 


i 

71 


0 1 0 \ 
-10 1 , 
0 - 10 / 


(80) 


and 


Jz = 


Based on these matrices and taking the parameter B as 
the one under estimation, one can easily check that the 
Hamiltonian H% satisfies the following equation 

[(F 3 x ) 2 - (x 2 + 4B 2 )]H*d e H 3 = 0, (82) 

which implies that if we choose fl = i/x 2 + IF 2 , the 
operator V defined in Eq. will be a conserved quantity. 
Specifically, V can be written in the form 

V = (x 2 - tt 2 ) Jz - 2 B X 1, (83) 

where X = |02)(20| + |20)(02|. Furthermore, the charac¬ 
teristic function Bb can be written as 

Ub = (X 2 / - t) J z - 2 B X 1 + ^ sin 2 (T) J Z X. (84) 


For the long-time limit, it reduces to 

—AB 2 t 

Hb = “ WXI ' (85) 

For the unitary parametrization processes with a purely 
initial state, the QFI for above Bb reads 

Fb = (pT^) {a2Jz) + 4 B V<a 2 f>, ( 86 ) 


where we have used the equality {J Z ,X} = 0. Further¬ 
more, the QFI can be simplified into 


F b 


/ 4 B 2 t V 

\X 2 + 4-B 2 J 


4B 2 X 2 


(J!) 


4 B 2 t 
X 2 + 4 B 2 



4 F 2 x 2 (F) 2 . 


(87) 


During the calculation we used the equality (J 2 ) = (X 2 ). 
The maximum value of QFI above is attained when 
(J 2 ) reaches its maximum value and (J z ), {X) vanish 
simultaneously. Denoting the initial state as |^) = 
ci| 02 ) + C 2 e l(f)2 \ll) + cse t(f>3 \20) with £ 1 , 2,3 real numbers, 
there are (J 2 ) = (X 2 ) = c 2 + c 2 , (J z ) = c 2 — c 2 and 
(X) = 2 C 1 C 3 cos 03. Utilizing these expressions, it can be 
checked that when ci = C 3 = l/y/2 and ^3 = 7 r/ 2 , all the 
conditions can be satisfied simultaneously, which implies 
that one optimal initial state here is a NOON-type state, 
i.e., 


l^o P t) = 4|(|02)+i|20». (88) 

The corresponding maximum QFI is 


F b , 


max 


= 4B 2 


4 B 2 t 2 . 

(x 2 +4 B 2 ) 2 X 


(89) 
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Single-mode cavity 



Figure 5: The schematic of a cavity optomechanics system. 
The left device is a single-mode cavity and the right one is a 
movable mirror. 


This expression shows that the maximum QFI will be 
square-enhanced with the passage of time. For the cou¬ 
pling constant x, Fb ,max does not change monotonously. 
For a small x> Fb, max increases sharply with the de¬ 
crease of x and the trend is totally reverse for a large 
X- The minimum value of Fb, max is attained around 
X 2 = 4(T>£) 2 / 3 — 4 B 2 . Thus, the value of x should be 
tuned carefully to avoid this regime during the measure 
of the transverse field. 


B. Canonical category 

The most obvious property of Hamiltonians in canoni¬ 
cal category is that they can be partly or entirely rewrit¬ 
ten via the canonical variables: x, p and the number 
operator. A typical case of canonical category is the op¬ 
tomechanical systems, which has been widely discussed 
as a novel artificial device m- a simple model for op¬ 
tomechanical systems is a single-mode cavity coupling 
with a movable mirror, of which the schematic is shown 
in Fig. [5] The total Hamiltonian can be written as [36] Gm 

H = uj a a^a + < jj^b — ga^a (b + , (90) 

where a, a\ b , b t are the annihilation and creation oper¬ 
ators for the cavity and mirror, respectively. cj a , uJb are 
the corresponding frequencies. The coupling strength g 
has the form 


9 = 


Ua 

l 


/ 


1 

2 rnuji) ’ 


(91) 


with /, m the length of the cavity and mass of the mirror. 
With the introduction of the quadratic operators x a (pp 
p a (b) defined as 


x a {p, T , Pa y^ 2 ‘ ^ 5 (92) 

Xb = 71 + ’ Pb = ~7Ti ~ b ^ ’ ( - 93 ' ) 


and the number operators N a = a^a, Nb = tfb, Hamilto¬ 
nian (90) can be rewritten into the form 

H = io a N a + io h N h - gN a x h , (94) 


where g = y/2g. In this case we take m or / as the 
parameter under estimation. Utilizing the commutation 
relations [N,x] = — ip , [N,p\ = ix and [x,p] = i, one can 
see that 


# x <9 m( z)# = iuj h g'N a p h , (95) 

(H x ) 2 d m(l) H = -ulg'NaXh+oJhgg'Nl (96) 

where we denote g' := d m ^g. From these equations, 

one can check that if we choose U = the operator 
V defined in Eq. is a conserved quantity. Its specific 
expression is 

V = t o b gg'N 2 . (97) 

As a matter of fact, the photon number in the cavity, 
i.e., N a is a conserve quantity, thus, it is natural that 
any exponentiation of N a is also a conserved quantity. 
Based on above information, the characteristic function 
reads 


^m(l) = g' N a [(t - u b f) x b + LO b (d t f)p b + gu> b fN a ] . 

(98) 

In the Fock space of the cavity, N a in above equation can 
be replaced by the average photon number n a , which is 
a constant, therefore, the characteristic operator can be 
simplified into 

Tl Q f 

{sin (u b t) x b + [1 — cos (u b t)] p b } . (99) 

LUb 

For a purely initial state, the QFI is in the form 

F m (i) = ^Fj{sin 2 (uj b t)(A 2 x b ) + [l - cos (uj b t)} 2 (A 2 p b ) 
+2 sin (ubt) [1 — cos (c^t)] cov ( Xb,Pb )} • (100) 

If the movable mirror is initially in the vacuum state, 
above expression of QFI reduces to 

F m(l) =(^fj [l - COS (C0 b t)}. (101) 


When t = {2k J rl)i:/ujb with k = 0,1, 2,..., the QFI above 
reaches its maximum value with respect to time, which 
is 4(1),max = (n a g'/ub) 2 - Contrarily, when t = 2kw/uj b , 
the QFI vanishes, the parameter cannot be estimated via 
Cramer-Rao inequality. This fact shows that the measure 
should not be performed at these time points. Moreover, 
the increase of photon number in the cavity can squarely 
benefit the estimation of l and m. 

For the mass m and the length /, the specific expres¬ 
sions of maximum QFI are 


n l u l 

4ra 3 / 2 co^ 


■ Z,max — 


n a^a 

ml 4 Luf 


( 102 ) 
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In both expressions above, tuning down the frequency ujq 
will help to improve the precision of l and m. Especially 
for the estimation of mass m, the decrease of uo 5 will show 
a dramatic enhancement of the precision. 

There are several other systems in this category, in¬ 
cluding a quantum harmonic oscillator in a classical field. 
The corresponding Hamiltonian is H = a + gaJ + 
g*a. In this case, the characteristic operator H is also the 
linear combination of operators x and p. The su(2) and 
canonical categories discussed above are representative. 
However, there are still systems out of these two cate¬ 
gories in which V is a conserved quantity. For instance, 
a two-level atom in a single-mode cavity with the Hamil¬ 
tonian H = u: a a^a + \wq(t z — g (a + <j z can also fit in 
the scenario discussed in this paper. 


V. ALTERNATIVE FORM OF QFI 


integrating form of T± can also be rewritten into an ex¬ 
panded form 


r±(G,6>)-f;g^(G x) n d e G. 

n =0 ' '' 


(108) 


For a thermal state expressed in Eq. ( 21 ), i.e., G = 
—/3H — In Z, where V is a conserved quantity, T± can 
be expressed by 

r± = fiV + (/ 1 ft 2 - p)d e H T 2 idpfi) H x deH. (109) 


where /j is defined as fi = 2Q,~ 3 [^Q — sinh(f O)]. Sim- 
ilarly with H in unitary parametrization process, T± is 
also the linear combination of V, the partial derive 3qH 
and its commutation between the Hamiltonian H . Thus, 
Iq for thermal states of all systems discussed previously 
can be calculated analytically. 


The classical Fisher information has more than one ex¬ 
tensions in quantum mechanics. Besides the traditional 
one discussed above, an alternative definition of quantum 
Fisher information is [28] 


Ie = 4Tr (de^fp) 2 . (103) 

For the unitary parametrization p(6) = U(6)poU^(6), 
this alternative form of QFI can be expressed by 


Ie = 8 Tr 


U 2 po - C^^/po) 2 , 


(104) 


where H is the corresponding characteristic operator. 
Similarly with the traditional expression, above formula 
is also determined by T~L and the initial state po • Recalling 
the spectral decomposition of po as po = Y^iL\ 
with M the dimension of the support of po, Eq. (104) can 
be written into 



(105) 

with (A 2 l~L)i the variance of T-L on the ith eigenstate of 
po- Similarly with the traditional form of QFI [24H27], 
Iq is also determined by the support of po- For a purely 
initial state, Iq = 8(A 2 H). 

For a general exponential form state po = exp (Gq), Iq 
is actually a correlation function, namely, 

I e = {T+{G,0)T-{G,0)), (106) 

where 

r ±(G,9):= [ e^^dgGds. (107) 

Jo 

When ei sGX doG is a real operator, T + = T_. Iq then 
reduces to (Tj_). Similarly with the SLD operator, above 


VI. CONCLUSION 

In summary, we discuss a general scenario in which the 
Hermintian operator V (defined in Eq. ([ 9 ])) is a conserved 
quantity. For the unitary parametrization processes, we 
provide analytical expression of the characteristic opera¬ 
tor which is totally determined by the Hamiltonian, 
the commutation between the Hamiltonian and its par¬ 
tial derivative, and the conserved quantity V. With the 
expression of H, we further give the expression of the 
QFI, calculate its maximum value and the correspond¬ 
ing optimal initial states. For the parametrized thermal 
states in this scenario, the SLD is the linear combination 
of V and 3qH. 

The scenario in this paper includes many specific phys¬ 
ical systems. As the application, we mainly focus on two 
categories: su(2) category and canonical category. In the 
su(2) category, we detailedly discuss the QFI in the ferro¬ 
magnetic two-spin system, the anisotropic two-spin XY 
model and a spin-one model. The characteristic oper¬ 
ator T-L in these systems can basically be expressed via 
su(2) generators. With the expressions of QFI, we locate 
the optimal initial states in these systems to access the 
maximum QFI. Meanwhile, the QFI for the parametrized 
thermal states of two-spin systems are also discussed. In 
the canonical category, we provide the QFI for a cavity 
optomechanical system. Increasing the photon number in 
the cavity or tuning down the frequency of the movable 
mirror will enhance the QFI. 

At the end of this paper, an alternative form of QFI 
is discussed in the scenario. Its formula for unitary 
parametrization processes is analytically given. For a 
general parametrized exponential state, we also provide 
the expression of this alternative QFI, which is a corre¬ 
lation function of T + and T_. In the scenario where V 
is a conserved quantity, T± is actually governed by the 
Hamiltonian, the parameter under estimation and the 
conserved quantity. We hope this paper could prompt 
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more and more researchers to study the connection be- No. 11475146. 
tween the conserved quantities and the quantum Fisher 
information, and search for various ways to enhance the 
parameter precision via conserved quantities. 
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